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Abstract. We discuss a clumpy model of obscuring dusty tori around AGN. Cloud- 
cloud collisions lead to an effective viscosity and a geometrically thick accretion disk, 
which has the required properties of a torus. Accretion in the combined gravitational 
potential of central black hole and stellar cluster generates free energy, which is dissi- 
pated in collisions, and maintains the thickness of the torus. A quantitative treatment 
for the torus in the prototypical Seyfert 2 nucleus of NGC 1068 together with a radiative 
transfer calculation for NIR re-emission from the torus is presented. 

1 Introduction: Dusty tori in the unified model of AGN 

Starting with the interpretation of spectropolarimetry of NGC 1068 by An- 
tonucci & Miller [1,2], the difference between type 1 and type 2 AGN is at- 
tributed to aspect-angle-dependent obscuration. The central continuum source 
is thought to be surrounded by a dusty molecular gas distribution, which is flat- 
tened by rotation into a torus or thick disk. The simplest unification scheme 
assumes that all Seyfert 2 nuclei harbor a Seyfert 1 core, so that the ratio of 
type Is to 2s measures the thickness of the torus. Present estimates range form 
1:4 [3] to 1:1 [4] when selected in the mid-infrared. The torus should therefore 
have a thickness HjR ^ 1. 

For more than a decade, clumpy models for the circumnuclear torus have been 
discussed [5]. It was recognized early on that dust will only survive in distinct 
clouds, which are embedded in the inter-cloud medium of the geometrically thick 
torus. This idea was nonetheless not included in the following radiative transfer 
calculations [6-9] of dust emission from these tori. Only recently, Nenkova et al. 
[10] developed a scheme which uses the dumpiness of the torus for an approxi- 
mative and statistical calculation of the thermal IR-emission from clouds, which 
are individually optically thick Ty > 40. This approach resolved some of the 
problems with earlier radiative transfer calculations. The toroidal distribution 
of atomic and molecular material is also inferred from observed Hi absorption 
and H2 emission, for example in the nucleus of NGC 4151 [11]. 

In two papers [12,13] we studied the consequences of the proposed dumpiness. 
It turns out that a scenario with distinct, quasi-stable clouds, which experience 
frequent cloud-cloud collisions, is unavoidable for geometrically thick torus. In 
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this paper we emphasis the importance of cloud coUisions and present the results 
of radiative transfer calculations for the torus in NGC 1068. 



2 Dusty clouds in the torus 

Cold, molecular and dusty clouds are the basic constituents of an accretion 
scenario of the torus. In spite of the complexity of their internal physics, we 
regard them as spherical clouds, which are sufficiently described by their radius 
rci and internal sound speed . The clouds must be self-gravitating to prevent 
dissolving in the inter cloud medium, which implies that the clouds free-fall 
time and the sound crossing time are approximately equal. This determines the 
volume filling factor 

of clouds, whore pQ — (j)vPc\ is the mean mass density in the torus and pci 
the density of individual clouds. The most important parameter of the torus 
in radiative transfer calculations is the vertical optical depth for intercepting a 
cloud 
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where /coU is the mean free path in the torus midplane (ZcoU is the local value) 
and ^ is the coefficient of the linear relation between po, pressure scale height 
H, and surface density S = ^poH , which follows from the vertical integration of 
the density. The radius rci of clouds can then be expressed as a function of r, 
Cs, and E. 

The mean free path is likely to be dominated by the largest clouds present. 
These clouds are more affected by tidal shear than smaller ones, which sets an 
upper limit to the possible cloud sizes. We expect that clouds accumulate at the 
shear limit when they experience increasing tidal forces while being accreted to 
the center. The upper limit is 

This is also a lower limit to the surface density 

Here M{R) is the total enclosed mass in stars and black hole at radius R and 
is the Keplerian circular velocity at that position. Both mass and optical depth 
will be dominated by the largest clouds and the right side of (4) will be a good 
measure of the true surface density. This implies a one to one correspondence 
between S and r. For the mass of individual clouds 
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Mcl < M{R)^^^ (5) 
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we find an upper limit, which is much smaller than the typical mass of molecular 
clouds in the ISM of normal galaxies. At a radius of 2pc with an enclosed mass 
of 10 ''Mq and an internal sound speed against gravitational cloud collapse of 
Cs = 2km/s gives a largest possible mass of 35 M0. Contrary to the cloud mass 
we find a lower limit for the hydrogen column density through an individual 
cloud 

iVH>»^ (6) 

which is about lO^^cm"^ in the environment close to the sublimation radius 
parameterized above. 

For the distribution of clouds in the torus we assume hydrostatic equilibrium 
for the vertical stratification. In [13] we used a modified isothermal distribution 
function of cloud velocities in an external potential, which includes a cut-off scale 
height Xh at which the density drops to zero. This leaves room for a wide polar 
funnel or an outflow cavity. We consider only cases, where the cut-off height is 
larger than the pressure scale height H. An example for the case of NGC 1068 
is shown in Fig. 1. The radial structure is derived from the stationary accretion 
scenario described below. 

3 Cloud collisions and accretion 

The unified model of AGN requires at least t ~ 1 for obscuration of the AGN for 
line of sights through the torus. As r is also a dimensionless collision frequency 

T ^ LOrJ fi, this implies that cloud-cloud collisions arc frequent in a torus. For 
anisotropic velocity dispersions of clouds Goldreich & Tremaine [14] derived an 
effective viscosity 

for angular momentum redistribution^. The required anisotropy can be deter- 
mined self-consistently for thin accretion disks and we use this limit also for the 
torus. 

Because in our scenario clouds in a torus are quasi-stable and only created 
or destroyed in cloud collisions, the momentum transfer from supernovae and 
stellar winds is inefficient for maintaining a large velocity dispersion of clouds. 
Accretion in the gravitational potential due to cloud collisions generate the free 
energy to balance collisional dissipation. 

Like in ordinary accretion disks, the effective viscosity from (7) allows mass to 
be accreted towards the black hole. The conservation law of angular momentum 
for the torus clouds with a viscous torque can be integrated once for stationary 
accretion 



^ The enhanced viscosity due to collective effects described by Wisdom & Tremaine 
[15] is unimportant for geometrically thick disks or tori with <C 1. 
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Fig. 1. Meridional cut through the probabihty density distribution of finding a dusty 
cloud in the torus. The distribution leaves room for an outflow along the polar axis 
and is derived from the model of Beckert & Duschl [13]. The spatial scale is in units 
of the dust sublimation radius. The mean number of clouds along a line of sight to the 
center drops below 1 for angels larger that 40° from the midplane {Z = 0). 

Special attention must be paid to the inner boundary at i?,; (index ?'). This 
boundary will be at the sublimation radius for diist (i? ^ 1 pc) , where neither 
torque nor shear will vanish. In addition the torque at the inner boundary is most 
likely not well described by the viscosity (7). Here M is the total mass accretion 
rate, which we assume to be constant at all radii throughout the torus. The outer 
boundary of our model is determined by the feeding process. The mass supply 
can be either a steady inflow driven by circummiclear starformation, discrete 
interactions with GMCs coming from large radii, or an inflow created by bars in 
the galactic potential. 

In cloud-cloud collisions a fraction ^(1 — e^) of the average relative kinetic 
energy of clouds is dissipated, where e is the coefficient of restitution, which 
approaches 1 for elastic collisions. We use e to parameterize cloud collisions and 
ignore the possibly more complex momentum redistribution in actual collisions. 
The energy dissipation appears as a sink term in the energy balance. For thin 
accretion disks with r ~ 1 we find e > 0.8, which is uncomfortably large. Smaller 
values of e are possible in geometrically thick tori, because energy advection can 
balance the increased dissipation in collisions. For a given set of parameters 
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(e, M, M{R)) the surprising result is tliat tori with high accretion rates are geo- 
metrically thick H/R ^ 2RM/ {cg M{R)) in spite of the smaller e of collisions. 
The elasticity, which leads to reasonably thick tori have e in the range 0.2-0.6, 
which might be provided by magnetic fields in the clouds. 

4 The Torus in NGC 1068 

For a detailed model of the torus one needs to know the mass distribution in 
the center. For the case of NGC 1068 wc collected the data in [13]. Some of the 
observed H20-maser spots [16] around the radio core component SI [17] trace 
a rotating disk or ring. The maser velocities appear inconsistent with Keplerian 
rotation around a central point mass. The rotation profile can be due to a massive 
disk [18], but we follow [19] and find a model with a black hole mass of 1.2 IQ^ Mq 
and a strongly concentrated stellar cluster with oc at large radii, a core 
radius of 0.32 pc and a core density of ^ = 5.25 lO^AfQ pc^'"*. Core radius and 
density are comparable to the nuclear stellar cluster in the center of our Galaxy 
[20]. The observational data from 0.3" to 4" from the nucleus are consistent with 
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Fig. 2. K-band brightness distribution of our radiative transfer calculations based on 
the method of Ncnkova ct al. [10] and the scenario of Bcckcrt & Duschl [13] for an 
inclination of 35° from the midplane. The spatial scale is in units of the dust sublimation 
radius. The contour scale of the surface brightness is logarithmic with a dynamic range 
of 2". 



this mass distribution. The central cluster extends out to 250 pc and the dusty 
torus lives in the gravitational potential of this cluster. 
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With the mass distribution, the model for individual clouds, and the cloud 

density distribution from Fig. 1 wc performed radiative transfer calculations 
based on the method of Nenkova et al. [10] . An example of the surface brightness 
distribution is shown in Fig. 2. The simulation demonstrate that the observed 
core component in Speckle images [21] in the NIR indeed measures the size of 
the sublimation radius to i? = 1 ± 0.4 pc. 

5 Conclusions 

We described a dynamical approach to the radial and vertical structure of a 
clumpy torus around an AGN following the ideas of Krolik & Begelman [5] and 
VoUmer et al. [12]. The vertical thickness is derived from an accretion scenario for 
a geometrically thick torus, which is based on cloud-cloud collisions. The inferred 
mass accretion rates for scale heights H > 0.5i? in the torus are larger than the 
Eddington rate for the central black hole and the model requires massive outflows 
from within the sublimation radius for dust. It is found that along typical lines 
of sight to the center (for angles < 45° w.r.t. the torus midplane) between one 
and ten clouds on average obscure the AGN. Only a small fraction of the mass 
accreted through the torus eventually reaches the black hole but is sufScient to 
generate the ionizing luminosity. An in depth analysis of cloud collisions in the 
environment of an AGN is required to test our model assumptions for e and the 
structure of individual clouds. 
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